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We have previously shown that spherically symmetric, inhomogeneous universe models can explain 
both the supernova data and the location of the first peak in the CMB spectrum without resorting 
to dark energy. In this work, we investigate whether it is possible to get an even better fit to the 
supernova data by allowing the observer to be positioned away from the origin in the spherically 
symmetric coordinate system. In such a scenario, the observer sees an anisotropic relation between 
redshifts and the luminosity distances of supernovae. The level of anisotropy allowed by the data 
will then constrain how far away from the origin the observer can be located, and possibly even 
allow for a better fit. Our analysis shows that the fit is indeed improved, but not by a significant 
amount. Furthermore, it shows that the supernova data do not place a rigorous constraint on how 
far off-center the observer can be located. 
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I. INTRODUCTION 



When analyzed within the framework of homogeneous 
FRW models, observations indicate that the universe is 
in a state of accelerated expansion. This is interpreted as 
proof that the energy content of the universe is presently 
dominated by some kind of dark energy with negative 
pressure Q,BI101[1I1I1I1]- However, for inhomoge- 
neous universe models this needs not be the case. A first 
step towards gaining an insight into how inhomogencitics 
affect the observables was made over two decades ago by 
Partovi and Mashhoon in Q and ^3- More recently, 
the interest in such models has grown increasingly as an 
alternative way to explain the appararent acceleration 
without introducing dark energy. A non-exhaustive list 
of references investi g ating this p ossibili ty is llll IT3I 
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In |l5| we showed that an inhomogeneous, but spheri- 
cally symmetric, matter-dominated universe model can 
easily explain the supernova Hubble diagram. Due to 
lack of a complete understanding of how perturbations 
behave in such inhomogeneous models, it is at present 
time not possible to calculate a theoretical CMB power 
spectrum to compare with that measured by the WMAP 
satellite [2^. However, one can still extract some infor- 
mation from the CMB, in the form of the location of the 
first peak in the power spectrum, that can be used to 
constrain the model. In [l^, it is shown that specific 
spherically symmetric models can account for both the 
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supernova Hubble diagram and the location of the first 
CMB peak. 

In such a scenario, the "accelerated expansion" is not 
a real effect, but is explained by the fact that the lo- 
cal expansion rate no longer is the same everywhere in 
space at a given time. Since the observations probe the 
expansion rate along the past light-cone, which at differ- 
ent redshifts corresponds to different points in space, the 
observed increasing expansion rate can just as well be 
explained by an appropriate spatially varying expansion 
rate. Qualitatively, the expansion rate would have to in- 
crease radially towards the observer. A specific model 
was found in Ref. [l^ where there is a very good agree- 
ment with both SNIa data and the location of the first 
CMB peak. In this model, the observer was located at 
the center of a large underdensity in an otherwise flat and 
matter dominated universe. Although the model is able 
to explain these data without the introduction of dark 
energy, it does so at the expense of the Copernican Prin- 
ciple, since the observer is located at a unique place in 
space. A more philosophically appealing scenario would 
be one where the observer does not necessarily have to 
be positioned exactly at the center of the underdensity. 

For an off-center observer in our model, the uni- 
verse appears to be anisotropic. In an earlier work 
|l6j . we showed that this anisotropy induces additional 
anisotropics in the temperature map on top of the in- 
trinsic anisotropies from the CMB. The farther away 
from the origin the observer is, the greater these induced 
anisotropies are. Thus, the observed CMB imposes a 
constraint on how far away from the origin the observer 
can be located. In we found the principal constraint 
to arise from the dipole, with a largest possible displace- 
ment from the origin of around 15 Mpc. Similarly, the lu- 
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minosity distance-redshift relationship will also become 
anisotropic. This might provide further constrains on 
the position of the observer inside the underdensity from 
an independent data set, namely the supernova observa- 
tions. 

In addition, it is in principle an interesting question 
whether the supernova data lend support to the idea of 
anisotropic surroundings, since the CMB indicates that 
the universe is extremely isotropic on large scales. This 
has previously been studied by Kolatt and Lahav 
and Bochner [3(|] . Although their approaches were some- 
what different, they both looked in essence for evidence 
of a lack of uniformity in the Hubble flow from a sta- 
tistical point of view. They both conclude that there 
is no statistically significant evidence for discarding the 
hypothesis of a homogeneous universe. However, the ad- 
vantage of our approach is that we can compare the data 
with a specific model, with actual predictions for the lu- 
minosity distance-redshift relation in different directions, 
rather than just treating angular variations in the super- 
nova data with statistical tools. Therefore, our method 
could yield interesting results despite the low number of 
supernova observations available. 

The structure of this paper is as follows. In Sect. HTl we 
present the solution to the field equations for the model in 
|l5| and look specifically at how the luminosity distance- 
redshift relation is altered for an off-center observer. In 
Sect. lIIII we calculate the theoretical luminosity distance- 
redshift relation and compare it with that inferred from 
the supernovae in the Ricss ct al. Gold Set. Finally, in 
Sect. Hvl we summarize our work and discuss the impli- 
cations of our results. 



II. THE LUMINOSITY DISTANCE IN AN 
INHOMOGENEOUS UNIVERSE 



The universe model we will focus on in this work is that 
described in ^5[. This is a spherically symmetric space- 
time, which is described by the Lemaitre-Tolman-Bondi 

(LTB) metric [3l|, |H H3] . The line element can be writ- 
ten as 



in terms of a conformal time rj 



l + /?(r) 



R 2 {r,t) [d6 2 + sin 2 Odcj) 2 } 



(1) 

where R{r, t) is a position-dependent scale factor, and 
/3(r) is related to the spatial curvature. For an exten- 
sive review of the properties of this model the interested 
reader is referred to Ref . j^il ■ 

For LTB models that contain only matter, the Einstein 
equations can be solved analytically. The exact form of 
the solutions depends on the sign of f3(r). For a model 
with a positive (3(r), corresponding to negative local spa- 
tial curvature, the solution can be written paramctrically 
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where a(r) is related to the density of matter. We must 
assume a(r) > for the matter density to remain positive 
everywhere. Furthermore, we have defined Rq = R{r, 0), 
interpreting t = as the time of last scattering, when 
photons decoupled from baryons. 

The scenario considered in [l5j was one where the ob- 
server is at the center of an underdense bubble. The 
angular diameter distance is then isotropic and is given 
by 



d A {z) = R(f) . 



(4) 



where f{z) describes the past light-cone of the observer at 
t = to. As indicated in the introduction, if the observer 
is placed at an off-center location, the distance measures 
will become anisotropic. The explicit effect this has on 
the expression for the angular diameter distance has been 
analyzed previously, first by Ellis et al. [35| , and also later 
by Humphreys et al. |3(| and Biswas et al. [23j|. We use 
the expression presented in |36| , which reads 



d\ sin 2 7 = .9 77 <? € £ 
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where g^ v is the observer frame, i.e. a frame centered on 
the observed in which the angular part of the line element 



dVL 2 = d 7 2 + sinV£ 2 



(0) 



The angles 7 and £ correspond to the polar and azimuth 
angles in this frame. 

In our case, such a frame can be constructed using the 
light-cones specified by (t,io,7,£), where t is cosmic time, 
and to is the time when the photons hit the observer at 
angles 7 and £. In order to simplify the calculations, but 
still keep the results completely general, we choose the 
z-axis in the direction of the off-center observer. The 
spatial coordinates of the observer in the reference frame 
defined by the metric Eq. (JJJ are then r = r b s and 6 = 0, 
with <f> being degenerate. Transforming the coordinates 
for the photon trajectories back to this frame from the 
observer frame, we get 

t = t (7) 
r = f(Mo,7) (8) 

(9) 
(10) 
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where the f and 9 functions are solutions to the geodesic 
equation as a function of t and the initial conditions rn, 
to and 7. Note that due to the axial symmetry about 
the z-axis, these functions do not depend on the azimuth 
angle £. A more thorough discussion of these solutions 
and their implications can be found in Ref . [Tfij . 

The metric g^ u in the observers local coordinate system 
is given by a simple coordinate transformation, i.e. 



which yields the following components 

2 
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Substituting these into Eq. © , we arrive at the following 
expression for the angular diameter distance 

(R') 2 _ fdr_\ 2 + fot 21 
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As a simple consistency check, we should recover the 
usual expression in Eq. Q in the limit where the ob- 
server is at the origin of the coordinate system. In this 
limit, the two angles and 7 coincide and the radial 
coordinate becomes independent of 7. This means that 
the two partial derivatives in Eq. (|15|) become and 1, 
respectively, and we end up with Eq. as expected. 

The distance measure probed by the supernova data is 
not the angular diameter distance, but rather the lumi- 
nosity distance d^. However, there exists a general re- 
lation [22} which allows us to relate these two distance 
measures: 



d L = (l + zfd A 



(16) 



Roughly, the two factors of (1 + z) take into account the 
reduction of energy per photon due to rcdshifting and 
the reduced arrival rate of incoming photons due to time 
dilation. 

Replicating the approach in 0], we solve the geodesic 
equations which determine the path of infalling photons 
relative to an off-center observer. Next, we use Eqs. I|15|) 
and l|16|) to obtain a theoretical prediction for the lumi- 
nosity distance-redshift relation in this scenario. How- 
ever, what we measure when we observe the supernovae 
is not the luminosity distance directly, but rather the ap- 
parent magnitude. These can be related to each other 
through the expression 



/i = 5 log! 
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where fi is the distance modulus, which is simply the 
apparent magnitude minus an absolute magnitude [3^ | . 
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FIG. 1: The relative matter density today as a function of 
physical distance from the center. 



III. RESULTS 



In the scenario which we considered in [l5j , the functions 
a(r) and /3(r) were parametrized as 



a(r) = H$r 3 
(3{r) = H 2 r 2 
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(18) 
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This corresponds to a smooth interpolation between two 
homogeneous regions where the inner region has a lower 
matter density than the outer region, thus describing a 
spherical bubble in an otherwise homogeneous universe. 
The parameter Hn = 100 /i ou tkms _1 Mpc _1 is the Hubble 
constant of the outer homogeneous region today, while an 
and Po are the relative densities of matter and curvature 
in this region. Furthermore, Aa and A/3 determine the 
difference in matter density and curvature between the 
regions, while r and Ar specify the position and width of 
the transition. In the original analysis we restricted the 
parameter space by imposing the constraints flo = 1 — ar, 
and A/3 = — Aa. We will keep these constraints in the 
present analysis too. 

Assuming that the observer was positioned at the center 
of the bubble, we then found a model that gave a good 
agreement with the Hubble diagram of observed SNIa 
and the position of the first CMB peak. The properties 
of this model are listed in tabled and the corresponding 
density profile of the undcrdensity is plotted in figure ^ 
The generalized matter density used in this plot is defined 



As explained in the preceding sections, we will investigate 
what effect moving the observer away from the center 
of the inhomogencity has upon the luminosity distance- 
redshift relation. To get an idea of how much the lumi- 
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Description 


Symbol Value 


Density contrast parameter 


An 


0.90 


Physical distance out to transition point [Gpc] 


D 


1.34 


Transition width 


Ar/r 


0.40 


Fit to supcrnovae 


2 

XSN 


176.2 


Position of first CMB peak 


s 


1.006 


Age of the universe [Gyr] 


to 


12.8 


Relative density inside underdensity 


ttrn,in 


0.20 


Relative density outside underdensity 




1.00 


Hubble parameter inside underdensity 


h in 


0.65 


Hubble parameter outside underdensity 


h out 


0.51 



TABLE I: The parameters and properties of the inhomoge- 
neous model which was studied in jlal . where the observer is 
located at the center of the inhomogeneity. 




0.5 



1.5 



nosity distance then varies across the sky, we have plotted 
in Fig. |2] the distance modulus as seen by an observer lo- 
cated at a physical distance D f, s = 200 Mpc from the 
center. The blue line in the plot represents the average 
distance modulus, while the shaded area is the RMS de- 
viation from the mean. Also plotted in the same figure 
are binned data points from the Riess et al. Gold Sample 
(RGS) of supernovae || and the best-fit ACDM model. 
It is evident that our model allow for relatively large 
variations in the distance modulus for rcdshifts z < 1 for 
off-center observers. One can therefore get a potentially 
even better fit to the data by placing the observer away 
from the origin, and taking into account the directions in 
which the various supernovae have been observed. 

The coordinates of all observed supernovae can be found 
at the CBAT website |4jj. The angular coordinates of 
the supernovae in the RGS are plotted in Fig. [3] The 
marker and color coding in the figure is as follows. Red 
circles correspond to supernovae with z < 0.5, blue pluses 
to 0.5 < z < 1.0, green stars to 1.0 < z < 1.5 and black 
crosses to z > 1.5. As we can see, the supernovae appear 
to be distributed relatively evenly across the sky except 
for in the proximity of the galactic plane. 

By moving the observer away from the center of the in- 
homogeneity, we add three additional degrees of freedom 
to the original model. These are the radial displacement 
from the center and the two angles that specify the di- 
rection of the displacement. We can treat these extra 
degrees of freedom in two different ways according to 
how we deal with the angles. First, we can average over 
all possible orientations of the underdensity. Basically, 
what we do then is throw away all information about 
angular dependence. If we didn't have any information 
about the direction in which the supcrnovae occurred, 
this would be the appropriate way to calculate the fit. 
The resulting \ 2 value as a function of the observer's po- 
sition is plotted as red circles in Fig. 0] It is clear that 
the fit becomes increasingly worse the farther away from 
the center the observer is placed. Thus, if one takes into 
account only the rcdshift and neglect information about 
the exact direct direction on the sky of the supcrnovae, 



FIG. 2: The average distance modulus for an off-center ob- 
server in our model. The shaded area represents the RMS 
deviation from the average across the sky. The data points 
and error bars are binned data from the Riess et al. Gold Set, 
while the red dashed line is the corresponding best-fit ACDM 
model. 




FIG. 3: The directions on the sky for 157 supernovae of the 
Riess et al. Gold Sample. Supernovae with z < 0.5 are 
marked with red circles, while supernovae with 0.5 < z < 1.0 
are marked with blue pluses, 1.0 < z < 1.5 with green stars 
and z > 1.5 with black crosses. 



the data disfavor an off-center observer. 

The other approach is to make use of the explicit angular 
dependence of the observed supernovae. Since this infor- 
mation is actually available to us, this would seem to be 
the more fitting way to do the analysis. The approach 
then is to minimize the y 2 with respect to the two angles 
for each radius. These minimal \ 2 values are plotted as 
blue crosses in Fig.^J It is evident from this plot that one 
can indeed reduce the residuals by moving the observer 
away from the center, if one also takes into account the 
specific directions in the sky of the supcrnovae. 

The minimized x 2 value is smallest for an observer lo- 
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FIG. 4: The x 2 fit to observed supernovae as a function of the 
observer's position. The red circles show the angle- averaged 
values, while the blue crosses show the angle-optimized values. 



cated at a physical distance of around D b s = 94Mpc 
from the center of the inhomogcncity. In Fig. we have 
plotted the x 2 seen by an observer at this distance as a 
function of the direction towards the center of inhomo- 
geneity . The x 2 reaches a minimum at Xmin = 174.9, 
with a direction towards the center of the inhomogeneity 
of (l,b) = (271°, 21°). 

As we discuss in |Tfij . such an off-center placement of 
the observer will necessarily induce an additional dipole 
in the temperature of the cosmic microwave background. 
CMB photons arriving at the observer's from the direc- 
tion of the center of the inhomogcncity will have traveled 
through a larger region with a high Hubble parameter 
compared to those from the opposite direction. They 
will therefore be more rcdshiftcd, and hence, appear to 
have lower temperatures. The COBE satellite [33 shows 
that the measured dipole in the background temperature 
points in the direction (l,b) = (264°, 48°). This means 
that the temperature is perceived to be higher in this di- 
rection compared to the opposite direction. Thus it ap- 
pears that the measured and the induced dipole for the 
best-fit off-center observer point in almost the opposite 
directions. However, this does not necessarily represent 
a problem, since the effect can be countered by an appro- 
priate peculiar velocity of our galaxy. We should there- 
fore not expect a correlation between the direction of the 
measured dipole and the direction towards the center of 
the inhomogcncity 



IV. DISCUSSION 



In this work we investigated the inhomogeneous but 
spherically symmetric model which we have explored prc- 



FIG. 5: The \ 2 residual for supernovae seen by an observer 
at a radial distance of D (, s = 94Mpc. The plot shows the 
residuals as a function of the direction towards the center of 
the underdensity. 



viously in and I n ^is model the observer is 

located inside an underdensity in an otherwise homo- 
geneous and flat universe, described by the Lemaitrc- 
Tolman-Bondi space-time. In |l5j we showed that it is 
possible for such models to explain both the supernova 
data and the position of the first peak in the CMB power 
spectrum. The aim of the current work was to learn 
whether it is possible to improve the fit to the super- 
nova data by moving the observer away from the cen- 
ter of the inhomogeneity, and also how much these data 
constrain such a movement. In addition to possibly pro- 
viding a better fit to the data, allowing the observer to 
be located at a different position than the dead center of 
the inhomogcncity is also more appealing from a purely 
philosophical perspective. After all, if we imagine plac- 
ing several imaginary observer randomly inside the un- 
derdensity, the chance of any of them ending up at the 
exact center vanishes. On the other hand, there is a fi- 
nite probability that an observer will end up inside any 
region with a finite radius. 

The universe no longer appears to be isotropic for such an 
off-center observer. This induced anisotropy is an effect 
that should be detectable in the observations. Specifi- 
cally for the supernova observations, this manifests itself 
in the form of an anisotropic relation between the lumi- 
nosity distance and the redshift. For such an observer, 
supernovae at the same redshifts but at different direc- 
tion in the sky, do not necessarily shine with the same 
brightness. 

We placed the observer at different radii from the center 
of the underdensity, and for each radius chose the direc- 
tion from the observer to the center to be that which 
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optimized the fit. In this way we were able to improve 
the fit compared to the isotropic observer. The fit was 
found to be optimal for a radius of D \, s = 94Mpc in 
the direction (l,b) = (271°, 21°) in galactic coordinates. 
However, the improvement turns out to be only slight. 
The minimal \ 2 f° r an off-center observer is 174.9 com- 
pared to 176.2 for an observer at the center. Although 
the x 2 has been reduced, the fit can not be said to have 
improved, considering that the off-center placement adds 
additional degrees of freedom. In fact, the x 2 per degrees 
of freedom is higher for all off-center observers compared 
to observers at the center. 

Looking at Fig. 0] we see that the % 2 is lower for off- 
center observers for radial distances out to about D = 
225 Mpc. From this we can conclude that anisotropies in 
the supernova data do not constrain very well how far 
away from the center the observer can be located. In 0| 
we looked at such constraints arising from the dipole in 
the CMB temperature. We found that the observer had 
to remain within a radius of around 15 Mpc for the dipole 
to be be in agreement with COBE measurements. It is 
clear that the current supernova data cannot improve this 
constraint any further. 



The main conclusion we can draw from our analysis is 
that the current supernova data do not offer any sub- 
stantial evidence for an off-center observer. Nor does it 
constrain very well how far off-center such observers can 
be located. This is partly due to the fact that there are 
too few supernovae in the sample. In the future, sub- 
stantially larger and better supernova samples will be 
available, such as those which will be provided by the 
Supernova Legacy Survey |4l|, and hopefully the Super- 
nova Acceleration Probe (SNAP) 01 ■ These might allow 
us to draw stronger conclusions regarding the anisotropy 
of the local universe. 
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